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Abstract 

Bond-dilution effects on spin-f/2 spin-gapped Heisenberg antiferromagnets of coupled alternating chains on a 
square lattice are investigated by means of the quantum Monte Carlo method. It is found that, in contrast with 
the site-diluted system having an infinitesimal critical concentration, the bond-diluted system has a finite critical 
concentration of diluted bonds, x c , above which the system is in an antiferromagnetic (AF) long-range ordered 
phase. In the disordered phase below x c , plausibly in the concentration region significantly less than x c , the system 
has a spin gap due to singlet pairs of induced magnetic moments reformed by the AF interactions through the 
two-dimensional shortest paths. 



1. Introduction 

Impurity effects on quantum Heisenberg antifer- 
romagnets with a spin gap due to the topological 
structure of the system have attracted much inter- 
est in relation to the impurity-induced antiferro- 
magnetic (AF) long-range order (LRO) observed 
experimentally [1,2,3,4]. There are two types of dis- 
orders, i.e., site dilution and bond dilution. Exten- 
sive theoretical and numerical works have inves- 
tigated properties of diluted Heisenberg antiferro- 
magnets in the one-dimensional (ID) [5,6] and two- 
dimensional (2D) systems. In the latter, however, 
the interchain interactions have been treated by 
the mean- field-type approximations [7] , so that the 
difference of these two types of disorders cannot be 
naturally introduced. In the present paper we have 
investigated bond-dilution effects on the spin- 1/2 
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spin-gapped AF Heisenberg system by treating the 
interchain interactions on the equal footing as the 
intrachain interactions by means of the quantum 
Monte Carlo (QMC) method. The result obtained 
clearly reveals the different aspects between mag- 
netic properties of the bond- and site-diluted sys- 
tems. 

The site-dilution effects have been already ex- 
tensively investigated [8,9,10]. If spins are ran- 
domly removed from the lattice, magnetic mo- 
ments are induced at sites neighboring the re- 
moved sites. We have called the domain of in- 
duced magnetic moments the 'effective spin'. They 
situate centered at sites next to each removed 
spin, and their extent is given by the correla- 
tion lengths, £,p' y , of the non-diluted system. The 
effective exchange coupling between two effec- 
tive spins centered at sites m and n is given by 

J mn oc (-l)l^-^+ilexp[-Z/(^^) 1/2 ], where 
I = \r m — r„ \ is the distance between the effec- 
tive spins [11,12]. There exist strong correlations 



between the effective spins retaining the stag- 
geredness with respect to the original lattice [10]. 
Therefore, the 2D site-diluted system has an AF 
LRO induced with an infinitesimal concentration 
of dilution at zero temperature. 

On the other hand, if a certain bond is removed 
from the lattice, the effective spins, peaks of which 
are located at both ends of the diluted bond, are 
induced. In this case, besides the same interaction 
Jmn between spins next to different removed bonds 
as in the site-diluted case, there exist AF interac- 
tions Jaf between two spins at both ends of a re- 
moved strong bond through the 2D shortest paths. 
The latter is estimated of the order of J' 2 , where 
J' is the strength of the interchain interaction and 
the strength of the stronger intrachain interaction 
is put unity. In the small concentration region with 
x <C l/£p£p> therefore it is expected that the sin- 
glet pair is reformed by Jaf- As x is increased, 
Jmn becomes dominant to Jaf- As a result, the 
AF LRO is expected to occur above a finite criti- 
cal concentration x c in the bond-diluted case. This 
scenario is clearly confirmed by the present work. 



2. Model and method 

We investigate the bond-diluted quantum AF 
Heisenbcrg model of coupled alternating chains on 
a square lattice described by the Hamiltonian 

H = X] e (2i,i)(2i+lj) S 2ij • S 2 i+lj 

+ a ^2 e (2i+l,j)(2i+2,j)^2i+lj ■ S 2 j+2,j (1) 

i,3 

+ J' X e (i,j)(i,j+l)^i,j ' SjJ+l , 

where 1 and a (> 0) are the AF intrachain al- 
ternating coupling constants, J' (> 0) the AF in- 
terchain coupling constant, and S^.j the quantum 
spin operator with S = 1/2 at site (i, j). Randomly 
quenched bond occupation factors {e(i,j)(fe,;)} in- 
dependently take either 1 or with probability 
1 — x and x, respectively, where x is the concen- 
tratin of diluted bonds. The pure system described 
by Eq. (1) with e(*,j)(fe,i) = 1 for all bonds is in 
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Fig. 1. Concentration dependences of the staggered mag- 
netization in the bond- and site-diluted [10] systems with 
a = J' = 0.5. There exists the critical concentration of 
x c ~ 0.05 in the bond-diluted system. The percolation 
thresholds on the bond and site processes arc denoted by 
x^, and Xp, respectively. All the lines are guides to eyes. 

a spin-gapped or an AF LRO phase depending on 
strengths of a and J' at zero temperature [13]. 
In the a = 0.5 system which we examine in the 
present work there exists the spin-gapped phase 
below J' c ~ 0.55. 

The QMC simulations with the continuous- 
imaginary-time loop algorithm [14,15,16] are car- 
ried out onLxL square lattices with the peri- 
odic boundary conditions. For each sample with 
a bond-diluted configuration, 10 3 ~ 10 4 Monte 
Carlo steps (MCS) are spent for measurement af- 
ter 500 ~ 10 3 MCS for thermalization. Sample av- 
erage for dilution is taken over 10 ~ 10 3 samples. 



3. Numerical results 

The x dependences of the staggered magnetiza- 
tion M s in the bond- and site-diluted [10] systems 
with a = J' = 0.5 are shown in Fig. 1. If the 
pure system with the spin-gapped state is diluted, 
the AF LRO is induced at a certain concentra- 
tion in both site- and bond-diluted cases. However, 
whereas the site-diluted system is driven into the 
AF LRO phase by an infinitesimal concentration, 
in the bond-diluted system there exists a critical 
concentration of x c ~ 0.05 even for a = J' = 0.5, 
which is the system near the critical point J' c ~ 
0.55 for a = 0.5 and x — 0. As the value of J' be- 
comes smaller, the region of the disordered phase, 
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Fig. 2. Real-space distribution of (a) the local static stag- 
gered structure factor and (b) the staggered correlation 
function in a fixed configuration with 15 diluted strong 
bonds for a = 0.5 and J' = 0.3 at T = 0.001. 

i.e., the value of x c increases due to the reduction 
of e.g., x c ~ 0.1 for J' = 0.3. 

For the coupled alternating chains considered in 
the present work there exist three types of bonds, 
whose strengths arc 1, a, and J', respectively. If 
a weaker bond with a or J' is removed, the sin- 
glet pairs on the strong bonds are not disturbed 
at all and so magnetic moments are not induced. 
However, if a strong bond is removed, magnetic 
moments are induced at both ends of the diluted 
bond. In Fig. 2(a) the real-space distribution of the 
local static staggered structure factor defined by 

S(i) = T.ji-^'^iSiSj) is shown for a = °- 5 
and J' = 0.3 on a 64 x 64 lattice from which 15 
strong bonds are removed. We can clearly see that 
magnetic moments are always induced in pairs at 
both ends of the removed bonds. This is a differ- 
ent situation from the site dilution for which iso- 
lated magnetic moments are induced next to each 
diluted site [10]. Furthermore, as shown in Fig. 
2(b), the real-space distribution of the staggered 
correlation function, C(i,0) = (-l)l r »" r ol {SfS§), 
exhibits no AF LRO, where (=(29,31)) is the 
peak site of an effective spin. (The reader may corn- 



Fig. 3. Spin gap proportional to J in the ex- 
tremely-low-concentration region for a = 0.5. A line with 
the slope of 1.5 is the guide to eyes. 

pare it with the AF LRO in the site-diluted case in 
ref. [10]). In such a small concentration region with 
x < where f£ = 3.452(1), Q = 1.8202(3) 

for a — 0.5 and J' — 0.3, the spins correlate only 
with the neighboring spins. In the concentration 
region of x ~ l/£p£p> on the other hand, the ef- 
fective spins practically overlap each other and the 
AF LRO is induced. 

Next we investigate whether or not there is 
the spin-gapped phase in the disordered phase by 
means of the moment method [16,17]. We have 
found that in the extremely-low-concentration re- 
gion there exists a spin gap, A, whose magnitude 
is independent of the value of x and proportional 
to J' 2 as shown in Fig. 3. This relation tells that, 
as mentioned in Introduction, two edge spins in- 
duced at both ends of a diluted bond reform a 
singlet pair due to the effective interaction Jaf 
(|Jaf| ~ J' 2 ) through the 2D shortest paths. We 
note here that since induced magnetic moments 
are located at sites of the different sublattices, the 
interaction is always AF. As J' 2 is increased, A 
deviates from the proportionality to J' 2 . This in- 
dicates that the localized singlet picture does not 
hold well for larger values of J'. 

In the inset of Fig. 4 the T dependence of the ef- 
fective Curie constant x( x )T for a = 0.5 and J' = 
0.3 is shown at x = 0.04, 0.02, 0.01, and from 
top. The uniform susceptibility x(0) exponentially 
decreases at low temperatures due to the spin gap 
with the value of 0.31415(2). In the main frame of 
Fig. 4, the differences (x(x) — x(0))T are shown. 



0.012 




Fig. 4. Temperature dependence of the effective Curie con- 
stant (x(x) - x(0))T for a = 0.5 and J' = 0.3. In the inset 
X(%)T are shown at x = 0.04, 0.02, 0.01, and from top. 
All the lines are guides to eyes. 

The plateau with the Curie constant proportional 
to x is clearly seen in these plots and (x(x)— x(0))T 
falls off exponentially at lower temperature than 
J' 2 = 0.09. This confirms again the existence of 
the spin gap due to Jaf at relatively small x. 



4. Summary and discussion 

When bonds are randomly removed in the spin- 
gapped phase, the AF effective interaction Jaf 
between two spins at both ends of a removed 
strong bond is induced in addition to the effec- 
tive interaction J mn between spins on ends of 
different removed bonds. Therefore, in such a low- 
concentration region that the strength of Jaf is 
larger than \J mn \ there exists the spin-gapped 
phase dominated by Jaf ■ 

The ground state of the diluted dimerized ID 
system with J' = and x > is in the quantum 
Griffiths phase, which is characterized by the finite 
correlation length and the gapless excitation. In 
the present 2D system with J' > our numerical 
results strongly suggest a gapless phase but with- 
out AF LRO just below x c . The problem will be 
discussed elsewhere. 
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